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Very-Near-Field  Plume  Model  of  a Hall  Thruster 

F.  Taccogna1,  S.  Longo1’2,  M.  Capitelli1,2 

1 Department  of  Chemistry,  University  of  Bari,  Italy 
2 IMIP-CNR  sect.  Bari,  Italy 

Plasmadynamic  phenomena  of  a stationary  plasma  thruster  (SPT-100)  plume  in  the  very-near-field  region  have  been 

studied  by  using  a two-dimensional  axysimmetric  numerical  code  based  on  a combination  of  Particle  in  Cell  (PIC/ 

+ ++ 

simulation  for  ion  component  (Xe  and  Xe  ) and  fluid  description  for  electrons.  The  model  can  reproduce 
qualitatively  several  experimental  observations. 
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1.  Introduction 

Stationary  plasma  thrusters  (also  called  closed 
electron  drift  thruster  of  the  Morozov  type)  [1]  are  ion 
thrusters  with  high  specific  impulse  and  low  thrust  that 
make  them  suitable  for  satellite  station  keeping  or  orbit 
transfer. 

A major  concern  in  the  use  of  these  devices  is  the 
possible  damage  their  plumes  may  cause  to  the  host 
spacecraft  and  to  communication  interference  of 
satellites.  Indeed  an  electric  thruster,  such  SPT-100,  in 
operation  produces,  besides  high  energy  ions  (Xenon 
for  his  inert  and  low  ionization  potential  property) 
responsible  of  the  thrust  and  electrons  emitted  to 
neutralize  the  positive  space  charge,  also  neutral 
propellant  atoms  and  low  energy  ions  created  by  charge 
exchange  (CEX)  collisions  between  ions  and  un-ionized 
propellant  (in  which  electrons  are  transferred)  or  by 
electron  impact  ionization  of  neutral  atoms.  They  are 
more  influenced  by  the  self-consistent  electric  fields 
that  cause  slow  ions  to  propagate  radially  and  to  flow 
upstream,  while  gaining  energy. 

Thus,  the  structure  of  the  plasma  plume  exhaust 
from  the  thruster  is  of  great  interest  and  several  models 
[2-5]  have  been  developed  simulating  the 
plasmadynamic  of  the  exhaust  plume  in  the  far-field 
region.  The  modelling  of  the  plasma  plume  of  electric 
thrusters  in  the  very-near-field  plume  region  (z<0.2  m 
from  the  thruster  exit  plane)  is  important  because  there 
are  many  new  aspects  to  take  in  account  as  plasma 
sheath  effect,  magnetic  field  effect,  non  isothermal 
electrons,  electron-neutral  collisions  and  so  on. 

2.  Numerical  model 

Due  to  particular  geometry  of  the  problem  each 
quantity  is  given  in  terms  of  axi symmetric  coordinates 
(r,z)  and  we  assume  no  variation  in  the  azimuthal 
direction. 

In  the  PIC  technique,  ions  are  considered  as  macro- 
particle and  their  charge  are  deposited  onto  a 
computational  grid  using  the  cylindrical  weighting 
functions.  Ions  are  loaded  into  the  simulation  at  each 
time  step  to  simulate  the  exit  flow.  The  ion  exit 
conditions  (radial  position  and  velocity  components)  are 
given  in  ref.  [5]  on  the  basis  of  fitted  experimental  data. 

In  order  to  know  the  electron  charge  density,  one 
differentiates  the  momentum  conservation  equation  and 
substituting  the  laplacian  of  the  electric  potential  with 
the  source  term  of  Poisson  equation  the  following 
differential  non  linear  equation  for  the  electron  charge 
density  results: 


(i) 

with  boundary  conditions  given  by  (see  Fig.  1): 

S2  ->  P,  = p~ 
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In  eq.  (1),  ve)t=vit  since  the  particle  flow  across  the 
magnetic  field  is  ambipolar  (current-less)  and  we 
consider  only  the  radial  component  of  the  magnetic 
field  (dipol  approximation). 


Fig.  1 - Scheme  of  the  simulation  domain. 

The  solution  is  obtained  by  the  iterative  Newton- 
Raphson  method  discretizing  the  laplacian  with  a 5- 
points  formula  in  cylindrical  coordinates  and  using  a 
successive-over-relaxation  (SOR)  scheme. 

Once  the  electron  charge  density  is  computed,  the 
electrostatic  potential  is  determined  by  solving 
Poisson’s  equation. 

Knowing  the  electron  density  we  can  solve  the 
electron  energy  conservation  equation  for  the 
temperature  where  we  neglect  the  unsteady  and 
convection  terms  and  we  consider  the  following 
ionization  energy  sink  term: 

a = KE;  (2) 

with  he  the  volumetric  ionization  rate  and  Ej 

the  Xenon  ionization  energy.  We  calculate  the  electron 
impact  ionization  rate  coefficient  k(Te)  on  the  basis  of 
e/Xe  cross-section  data  [6]  fitted  to  the  following 
polynomial  expression: 
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k(T_)=  4.5243  x 10'”  -4.4952  x 10"” Te  + 1.1776  x 10'”  fr  - 
255 1 5 x 1 0 -'T/  - 2 .3932  x 1 0-|K  7;.J  + 9 .7869  x I O'20  Tj 

(3) 

Substituting  the  ideal  gas  equation  we  have 

k,V27;,  + Vk(,-V7>-£.  (4) 

where  the  electron  thermal  conductivity  is  included  by 
the  following  fitting  formula  (including  the  electron 
collision  effect)  [7]: 


k„  = AT? 


with 


\A  = 6.9  x 1(T12 
1 A = 2.08  x 10~3 


w = 5/2, 
n = 1/2, 


Eq.  (4)  is 
conditions: 


solved  with  the  following 


5,,^  — > Te  = const 
dT 

or 

L2,Ly,LA  T = T 


l <2eV 
Tc  >2eV 

(5) 

boundary 


and  we  come  back  to  the  electron  momentum 
conservation  equation  until  we  reach  convergency. 

The  ions  are  moved  under  the  influence  of  the  self- 
consistent  electric  field  integrating  numerically  the 
equations  of  motion  using  the  leap-frog  algorithm  . 

Collisiona!  processes  are  modeled  between  move 
steps  using  a TPMC  method  [8]  to  simulate  Xe-Xe  and 
Xe-Xe  (momentum  and  charge  transfer)  collisions, 
and  a volumetric  production  method  to  simulate  neutral 
ionization  by  electron  impact. 

3.  Results 

The  behavior  of  the  plasma  potential  is  reported  in 
Fig.  2.  We  can  see  the  potential  wall  structure 
developing  spontaneously  that  capture  ions  to  neutralize 
electrons  trapped  by  magnetic  field.  In  fact  we  note  a 
peak  near  the  centerline  at  the  small  axial  position  that 
is  attributed  to  electrons  confined  by  magnetic  field 
cusp  formed  by  the  thruster  magnetic  coils. 


Figure  3 shows  the  radial  profiles  of  electron 
temperature  measured  [9]  at  different  axial  locations 
The  peak  structure  in  front  of  the  discharge  channel 
(28mm<r<50mm)  decreases  in  magnitude  and  broadens 
as  the  electrons  move  away  axially  from  the  thruster  in 
agreement  to  the  fact  that  electrons  continue  to  cool 
down.  We  note  a discrepancy  between  the  model  and 
experimental  values  in  the  very  low  radial  and  axial 
positions  due  to  neutralizer  position  effect.  The  high 
experimental  values  of  temperature  at  radial  distances 
greater  than  70  mm  for  z=10  mm  can  be  due  to  a 
reduction  of  thermal  electrical  conductivity  caused  by 
electron-ion  collisions. 


Radial  Posriion  (mmj 

Fig.  3 - Measured  [9]  and  computed  electron 
temperature  as  a function  of  radial  position  for  different 
axial  position. 

Acknowledgements 

This  paper  was  partially  founded  by  the  AS!  under 
contract  ARS  IC-1 42-00-0.  The  authors  thank  Gianpiero 
Colonna  (IMIP,  CNR)  for  providing  the  set  of  e/Xe 
cross-sections. 

References 

[1]  V.  V.  Zhurin  et  al,  Plasma  Source  Sci.  Technol.  8 
(1999)  R1-R20. 

[2]  D.  Y.  Oh  et  al,  J.  Propulsion  and  Power  15 
(1999)  345-357. 

[3]  I.  D.  Boyd,  D.  B.  VanGilder,  M.  Keidar,  J. 
Spacecr.  Rockets  37  (2000)  129-136. 

[4]  F.  Taccogna,  S.  Longo,  M.  Capitelli,  A1AA  paper 
2000-2345  (2000). 

[5]  F.  Taccogna,  S.  Longo,  M.  Capitelli,  J.  Spacecr. 
Rockets  39  (3)  (2002)  409-4 1 9. 

[6]  G.  Colonna,  L.  D.  Pietanza,  M.  Capitelli,  AIAA 
Paper  2000-2349  (2000). 

[7]  R.  Samanta  Roy,  Ph.  D.  thesis , MIT  (1995). 

[8]  K.  Nanbu,  Y.  Kitatani,  J.  Phys.  D:  Appl.  Phys. 
28(1995)324-330. 

[9]  S.  Kim,  J.  E.  Foster,  A.  D.  Gailimore,  AIAA 
paper  96-2972  (1996). 


Fig.  2 - Plasma  potential  surface. 
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